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Abstract. The Chamaeleon I dark cloud (Cha I) has been mapped in C ls O ( J = 1 — 0) with an angular resolution of 1' 
using the SEST telescope. The large scale structures previously observed with lower spatial resolution in the cloud turn into a 
network of clumpy filaments. The automatic Clumpfind routine developed by Williams et al. 1994 is used to identify individual 
clumps in a consistent way. Altogether 71 clumps were found and the total mass of these clumps is 230 M©. The dense 'cores' 
detected with the NANTEN telescope l Mizu no et al. 199 9 1 and the very cold cores detected in the ISOPHOT serendipity survey 
<T6th et aT~2 000 1 form parts of these filaments but decompose into numerous 'clumps'. The filaments are preferentially oriented 
at right angles to the large-scale magnetic field in the region. We discuss the cloud structure, the physical characteristics of the 
clumps and the distribution of young stars. The observed clump mass spectrum is compared with the predictions of the turbulent 
fragmentation model of Padoan & Nordlund 2002 An agreement is found if fragmentation has been driven by very large-scale 
hypersonic turbulence, and if it has had time to dissipate into modestly supersonic turbulence in the interclump gas by the 
present time. According to numerical simulations, large-scale turbulence should have resulted in filamentary structures as seen 
in Cha I. The well-oriented magnetic field does not, however, support this picture, but suggest magnetically steered large-scale 
collapse. The origin of filaments and clumps in Cha I is thus controversial. A possible solution is that the characterization of 
the driving turbulence fails and in fact different processes have been effective on small and large scales in this cloud. 
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1. Introduction near Ced 110 and Ced 112 is provided by observed molec- 
ular outflows JMattilaet al. 19891 IPrusti et al. 199U and two 

The dark cloud Cha I is one of the nearest low-mass star mm con tinuum sources pelgurth et al. 1996). Th e research on 

forming regions. The distance 150 pc to the cloud has been cha l until mQ is rev iewed bv lSctiwartz 19911 A short sum- 

g recently determined using the Hipparcos satellite data by ma ry of later large scale studies is given in the following. 
Knude & H0g 1998 The elongated cloud with the apparent 

dimensions ~ 0.5° by ~ 2° is easily recognized on the A search for new y° un S stellar ob J ects in Cha 1 has 
ESO/SRC sky survey plates. A visual extinction map based been conducted by |Cambresy et~ 998 using DENIS UK 
on star c ounts on the ESO/SRC sky survey plates is pre- photometric data. Most sources were distributed in a sim- 
sented in IToriseva & Mattila 19851 and a near-IR extinction ilar wa y t0 the known member candidates in the central 
map from DENIS UK star counts by |Cambresy et al. 1997| and northern parts of the cloud. Only few young star can- 
Several signposts of past and present star formation are ob- didates w ere found in the southe rn part. A deeper JHK 
served in the cloud. Three visual ly bright reflecti on nebulae stud y b Y |G6mez & Kenyon 2001| confirmed these results. 
(Ced 1 10, Ced 1 1 1 and Ced 1 12; Icidgrblad 1946 1 and an in- New candid ate member stars have be en found using L-band 
frared reflection nebula (the IRN) are associated with Cha I. photom etry <Kenyon& Gomez 2001} and JHK variability 

. . . studies (Carpenter et al. 2002 1. Using ISOCAM observations 

More than hundred pre-main-sequence stars have been — — . — Im—LLf-r- ■ — 1 ° . 

c j • sl. i j r. • j i , . t., a Persi et al. 2000 found clusters of sources with rmd-IR excess 

found in the cloud. Besides clusters around the visible reflec- K 1 

. i . ., around the three reflection nebulae. A number of sources with 

tion nebulae, these objects are scattered around its western mid . , . . , 

. . , „ . , „ . „ . mid-IR excess is seen also in the very eastern and western parts 

part and the northern part. Evidence tor ongoing star formation , , , , , • , , . , ,. „ , 

of the cloud but only one in the dense region southwest of Ced 

111. 

Send offprint requests to: L. Haikala 

* Based on observations collected at the European Southern The multicolour IRAS images l |Boulanger et al. 1990| 

Observatory, La Silla, Chile |Boulanger et al. 19 98 1 outline the dense parts of Cha I. FIR 
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ISOPHOT observations of the southern part are reported in 
Haikala et al 1998. Besides the emission from the strong IR 
point sources in the reflection nebulae Ced 110 and 111 and 
the IRN, extended FIR emission is seen to coincide with the 
extinction features mapped by Toriseva & Mattila 1985] and 
|Cambresy et al. 1997| In particular two extended dust emis- 
sion maxima are seen southwest of Ced 111. These dust emis- 
sion maxima were detected also by the ISOPHOT 170 /xm 
Serendipity Survey observations of Cha I (ITofh et al. 20 00 1. 
However, up to now no FIR point sources have been detected 
in these maxima. 

Compared to the major star formation regions in the north- 
ern sky, radio molecular line observations of the Chamaeleon 
region have been few. The first large-scale radio molecular 
line mappings of the cloud were made with the Parkes 64- 
m telescope in the 6 and 18 cm transitions of H 2 CO and 
OH, respectively ( Tori seva et al. 19851 . The 1.2-m Columbia 
Southern Millimeter- Wave Telescope at Cerro Tololo was used 
to map Cha I in the C 18 ( J = 1 — 0) (IToriseva et al. 1990i 
and the 12 CO ( J = 1 - 0) ( Boulange r et al 1998) . Due to the 
small size of the telescope, the spatial resolution was moder- 
ate (9' HPBW). First observations in the 1' resolution range 
were obtained when the SEST telescope started observations 
in 1988. As the beam of the telescope is 45" to 55" at 3 
mm, only small areas were mapped (Mat tila et al. 1 989 1. The 
Nagoya 4-m survey telescope, NANTEN, located at the Las 
Campanas Observatory, was used for mapping the Chamaeleon 
region in the 12 CO, 13 CO and C 18 Q (J = 1-0) transition s 
(Miz uno et al. 20011 IMizuno et al. 19991 IMizuno et al. 1998i . 
The beam size of the NANTEN telescope is 2.7 arcminutes but 
the observing grid was 8', 4' and 2', respectively. 

This paper reports the mapping of the dense parts of the 
Cha I cloud in the C 18 (J = 1 — 0) line with an angular res- 
olution of 1'. These observations enable more detailed studies 
of the cloud structure than has been possible by using optical 
or NIR star counts (statistical limitations) or by the smaller size 
radio telescopes (diffraction limited beam sizes of several arc 
minutes). The relatively high frequency resolution of 43 kHz 
(corresponding to ~ 0.12 kms -1 at 109 GHz) makes it pos- 
sible to extend the analysis also into the velocity space. An 
automated approach to data analysis, the "Clumpfind" routine 
developed by Willi ams et al. 19941 has been used to analyze the 
(x, y, v) data cube derived from the observations. 

The observations and the data reduction are described in 
Sect. 2. In Sect. 3 we describe the cloud structure and results 
of the clump finding procedure, and discuss the new results in 
the light of previous observations concerning the signposts of 
star formation and magnetic fields in the cloud. A comparison 
with the NANTEN cores and clumps identified in this paper is 
presented in Sect. 4, and notes of some individual regions of 
special interest are presented in Sect. 5. In Sect. 6 we discuss 
the clump stability and mass spectrum. Finally, in Sect. 7 we 
summarize our conclusions. 

2. Observations 

The observations were conducted during various observ- 
ing runs from 1988 to June 1996 using the Swedish-ESO- 



Submillimetre-Telescope SEST on La Silla observatory in 
Chile. A 3 mm dual polarization single sideband (SSB) 
Schottky receiver in the frequency switching mode (with a 
throw of 7 MHz) was used to observe the J = 1 — tran- 
sitions of C ls O, 13 CO and C 17 at 109, 110 and 111 GHz, 
respectively. The SEST high resolution, 2000 channel acousto- 
optical spectrometer with an 86 MHz bandwidth and a chan- 
nel width of 43 kHz was split into two bands to measure the 
both polarizations simultaneously. The last C ls O observations 
(1995,1996) were done using the SEST 3mm SSB SIS receiver. 

Calibration was achieved by the chopper wheel method. As 
the observations are spread over many observing runs with dif- 
ferent weather conditions the system temperatures varied ac- 
cordingly. Typical values for the effective SSB system temper- 
atures outside the atmosphere ranged from 350 K to 450 K (for 
the Schotky receiver) and from 200 K to 300 K (SIS receiver) 
for the C ls O measurements, and from 300 K to 350 K for the 
13 CO measurements. 

During the different runs the same reference position was 
always observed to monitor the calibration. Pointing was 
checked every 2-3 hours towards the nearby SiO maser source 
U Men. We estimate the pointing accuracy to have been bet- 
ter than 10" during the observing runs. The focusing was done 
using a strong SiO maser. 

C ls O Altogether 1836 positions in the cloud were ob- 
served using a map step size of 1' and integration time of 60 
seconds. The observations cover ^0.6 square degrees. The map 
is not fully sampled since the SEST beamsize at this frequency 
is 45". We believe however that the clump detection is mainly 
limited by the noise level reached, and not by the slight under- 
sampling. The 1' step size corresponds to ~ 0.044 pc at the 
distance of 150 pc to the cloud. The (0, 0) position was arbi- 
trarily chosen to coincide with star T39 near the cloud centre 
(J2000.0: ll h 09 m 13!3, -77°29'12"). The C ls O observations 
of the Ced 110 and 112 regions have already been reported in 
IMattila et al.T989l 

The median RMS noise of the C ls O spectra, after folding 
and baseline fit (typically second order baseline was used), is 
0. 1 K. Because of different observing conditions the noise goes 
occasionally up, and the maximum level is 0.22 K. For 90% of 
the spectra the RMS noise is smaller than 0.15 K. 

The integrated C 18 ( J = 1 — 0) intensity map shown in 
Fig- n serves as a finding chart for the identified clumps (see 
Section 3) and some prominent objects in the cloud. The lo- 
cations of the cores identified by Miz uno etal. 19991 and vis- 
ible reflection nebulae Ced 110, Ced 111, and Ced 112, and 
the IRN are indicated. Also shown are the locations of two 
Class sources (Cha-MMSl and Cha-MMS2) detected by 
|Reipurth et al. 1996| A velocity channel map over 20 channels 
of 0.12 kms" 1 in width is shown Fig. [2] The velocities are 
indicated in the top of each panel. 

13 CO The central and southern parts of Cha I were ob- 
served in 13 CO ( J = 1 — 0) . Observing was done in a similar 
manner as the C 18 observations but with a 2' map step size. 
The number of observed positions is 296 and they cover an area 
of ^0.25 square degrees. The integrated T^ map is shown in 
Fig. |3] The locations of Ced 1 10, Ced 1 1 1, and the IRN are also 
indicated. 
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Fig. 1. Integrated C ls O intensity map over the LSR velocity range 2.9 — 5.9km s -1 of the Chamaeleon I cloud. The lowest 
greyscale level is 0.2 K km s^ 1 and the stepsize is 0.3 K km s _1 . The lowest contour level is 0.8 K km s~ x . The crosses indicate 
the positions of Miz uno et al. 1999l C 18 Q cores 3,4,6 and 5 (from north to south). The locations of the reflection nebulae Ced 112, 
Ced 110 and Ced 111 (from north to south) are indicated with open squares. The locations of the Class protostar candidates 
Cha-MMSl and Cha-MMS2 are marked with filled circles. TRN' marks the location of the Infrared Reflection Nebula.The 
numbers identify the clumps found in the cloud (see Table. 
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Fig. 2. Maps of the C ls O intensity in each of the central channels in the Chamaeleon I cloud. The LSR velocity in each panel 
is indicated in the upper corner of the panel. The highest intensity in the figure is ~ 2.2 K 
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Fig. 3. The integrated 13 CO(J = 1 - 0) intensity map 
of the central part of Cha I. The lowest greyscale level is 
2.2 K km s _1 and the stepsize is 0.6 K km s _1 . The markers are 
the same as in Fig.[0 The box in the centre of the figure outlines 
the area displayed in Fig.|4] 



Fig. 4. Sample 13 CO and C ls O J = 1 — spectra observed in a 
6' by 2' region indicated in Fig. [3] The maximum 13 CO T^ in 
the figure is 3K. The offsets relative to the map centre position 
(R.A. = ll h Q9 m 13!3, Dec. = -77°29'12") in arc minutes are 
shown on the axes. 

C 17 Pointed C 17 ( J = 1 - 0) observations were 
made towards a few C ls O maxima. The integration time was 
20 minutes, and the frequency throw was set to 12 MHz due 
to the hyperfine structure of the line. Longer C ls O integrations 
were made towards the same positions. In Fig. [5] we present 
C ly O and C ls O lines observed in three positions in the cloud. 

3. The cloud structure 

The integrated C 18 and 13 CO ( J = 1 - 0) intensity maps 
presented in Fig. [T] and Fig. |5] show similar features on a 
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Fig. 5. A comparison of the C 18 ( J = 1 — 0) (histogram) 
andC 17 ( J = 1 — 0) , F = 5/2- 5/2 line (continuous line) 
intensities in selected C 18 maxima. The C 17 intensity has 
been multiplied by 9, which corresponds to the expected inten- 
sity ratio of the C 18 ( J = 1 - 0) line to this C ir O hyperfine 
component in the optically thin case. 



large scale. As the 13 CO mapping was done using 2' step 
size it can not show the cloud structure in smaller detail. 
A closer look reveals significant differences between C 18 
and 13 CO. For example, the extended C ls O maximum near 
ll h 03 m 30 s , -77°37' is not visible in the 13 CO map. Another 
example is given in Fig. 0] The intensities of the C ls O and 
13 CO ( J = 1 — 0) lines towards the shown positions do not 
correlate with each other. The reason for this failure is that the 
13 CO line becomes optically thick in the direction of a high 
column density. The case of 12 CO is even worse as it is clearly 
strongly self absorbed in the line centre almost everywhere (see 
eg. Figures 3 and 19 in Mani la et al. 19891 . 

The comparison of the rare CO isotopomer C 17 and C ls O 
lines towards three positions in the cloud is shown in Fig. [5] 
The two transitions are well correlated and their intensity ratio 
(when the three C 17 ( J = 1 — 0) components are added) is 
3.6, which corresponds to the abundance ratio in the local ISM 
(Wilson & Rood 1994i. This comparison shows that both lines 
are likely to be optically thin and that the C ls O( J =1-0) line 
can be used to trace the total CO column density distribution. 

A conspicuous feature in the C ls O maps is the presence of 
several clumpy filaments. The most prominent of these is the 
northern arm, a 40' long north-south oriented structure reach- 
ing down from Ced 1 12. In the centre and in the south, there are 
three almost parallel northeast-southwest oriented filaments, 
and an arc almost 40' in length, stretching from the neighbour- 
hood of the IRN towards the northern arm. When investigated 
in detail these larger entities break up into smaller clumps. 
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The notable northeast-southwest oriented filaments break into 
small, often north-south oriented, clumps. 

3.1. Clump identification 

We have analyzed the small scale structure seen in the C ls O 
map by using the automatic routine "Clumpfind" developed by 
IWilliams et al. 1994I which assumes that intensity maxima in 
the spectra correspond to localized density enhancements or 
'clumps'. The C ls O(J =1 — 0) spectral line map was inter- 
polated to a three-dimensional T^(x,y,v) data cube, where 
the pixel size was 30" in the position coordinates (one half 
of the map grid spacing), and one spectrometer channel width 
(0.12km s _1 ) in the radial velocity. The C ls O emission from 
the cloud is confined to the velocity range 2.7 — 6.1 kms -1 
and only the central 29 channels corresponding to this range 
were included to the data cube. The maximum dimensions of 
the map in R.A. and Dec. directions are 41' and 96', respec- 
tively, and the dimensions of the data cube are 82 x 192 x 29. 

According to the analysis of Willi ams et al. 199"4*l the inten- 
sity stepsize, AT, for the Clumpfind routine should be set to 
two times the maximum noise level. Instead of the absolute 
maximum, we used here the 'high' value, 0.15 K, below which 
the noise level remains for 90% of the spectra, and accordingly 
the intensity stepsize was set to 0.3 K. This was used also as 
the intensity threshold for the clump search and the total cloud 
mass estimate. 

With these parameters the Clumpfind routine identifies 7 1 
clumps in the dataset. The data cube ("Chal.fits") and the 
clump identification file ("Chal.fits. elf"), which is an output of 
the Clumpfind programme, will be made electronically avail- 
able as FITS files. An IDL routine called 'cLsurf.pro' written 
by us to help in inspecting the cloud structure will also be avail- 
able. This small programme shows clumps as isointensity sur- 
faces in the (x, y, v)-space using the IDL Object Graphics and 
the 'xObjView' interface for zooming and rotating these sur- 
faces. 

The physical characteristics of the clumps are presented in 
Tabled The columns of this Table are: (1) the clump iden- 
tification number; (2) and (3) the Galactic coordinates of the 
clump centre of mass; (4) and (5) R.A. and Dec. (2000.0); (6) 
and (7) the LSR velocity of the line peak and the FWHM of the 
velocity dispersion in the clump; (8) the clump half -power ra- 
dius; (9) the maximum C 18 column density; (10) the clump 
mass with an error estimate based on the spectral noise; and 
(11) the ratio of the gravitational potential to kinetic energy in 
the clump. The numbering corresponds to order in decreasing 
peak intensity. Clumps 1-6 are detected at the level 7 x AT 
(14cr), clumps 7-13 at the level 6 x AT, clumps 14-20 at the 
level 5 x AT, clumps 21-36 at the level 4 x AT, clumps 37-52 
at the level 3 x AT, and clumps 53-71 at the level 2 x AT (4cr). 

The division of the data cube into clumps in the central 
and northern parts of the cloud is demonstrated in Fig.|6] The 
clumps identified by the Clumpfind routine are plotted in the 
figure as surfaces in the (x, y, v) space, where x and y are the 
R.A. and Dec. offsets from the (0, 0) position, and v is the LSR 
velocity. The total velocity spread in Cha I is only 2.3 kms -1 



but sometimes two clumps can be separated along the same line 
of sight. 

The C 18 column densities (for each x,y,v pixel) were 
estimated by assuming optically thin emission, LTE with the 
excitation temperature T cx = 10 K, and a beam-source cou- 
pling efficiency, ijc, of 0.8 (half-way between 7/mb and ijMoon)- 
These were converted to mass column densities by assum- 
ing that the column density ratio N(C 18 0)/N(K 2 ) is 2 10~ 7 
(Wil son & Rood 19941 . and that the gas contains 10 % He. The 
masses were derived by summing up all pixels within a clump. 
The distance assumed in the mass calculation is 150 pc. For the 
gravitational potential energy estimates, it was assumed that the 
clump diameter in the radial direction (z) is equal to the smaller 
of the diameters in x and y directions. The mutual gravitational 
potential energies of all (x, y, z) pixel pairs were summed up 
and the internal potentials of pixels (small contribution) were 
added by approximating them with homogeneous spheres. The 
kinetic energies were also calculated pixel by pixel, and they 
include the contributions of systematic and turbulent motions 
and the internal thermal energy. 

The total mass of the cloud (as traced by C 18 above the 
threshold 0.30 K) is 230 M Q . The derived clump masses range 
from 0.5 to 12 M Q with the median 2.9 M©. The geometric 
mean radii lie in the range 0.08 — 0.21 pc. The ratios of the 
gravitational to kinetic energies are between 0. 1 and 1 .2, and 
the median ratio is 0.4. The masses and the stability of the 
clumps will be briefly discussed in Sect.|6] 

The largest uncertainty of the mass estimates is related to 
the the N (C 18 0) /iV(H2) conversion factor. The value we have 
used is close to those derived by |Frerking et al. 1982| in Taurus 
and[Harjunpa a & Mattila 1996 in Cha I, and thereby also con- 
sistent with the factor used by IMizuno et al. 1 999 The con- 
version factor is likely to change towards the centres of dense 
clumps due to CO depletion. However, as the column densi- 
ties are generally modest in this cloud and most dramatic de- 
pletion effects are localized to very high density regions (e.g. 
ICaselli et al. 19 99 1 we believe that CO depletion does not cause 
significant errors to the clump statistics. 

3.2. Star formation and cores 

The locations of known and likely members and candi- 
date pre-main-sequence stars of the Cha I Association are 
shown in Fig. [7] projected on the C ls O molecular line 
map. The stellar objects were selected from ISchwartz 19911 
(Optical candidate members, black triangles), P ersi et al. 2 000 
(known members with NIR excess (red squares), mem- 
ber candidates with NIR excess (open red triangles)), 
Cambresy et al. 1998 (new Young Stellar Object candidates 
(asterisks)) and Gomez & Kenyon 2001 (candidate pre-main- 
sequence stars (open circles)). Many of the objects listed in 
ISchwartz 199T1 were replaced by known members with NIR 
excess from P ersi et al. 20001 A number of investigations have 
been conducted in the direction of the three reflection nebulae 
known to be locations of active star formation and more candi- 
date members could be selected from them. We think, however, 
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Table 1. Clump properties. 
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that the objects indicated in Fig.0reflect well the general dis- 
tribution of young stellar objects in Cha I. 

Concentrations of young stars are seen in the direction of 
the three reflection nebulae, Ced 112 in the north, Ced 110 in 
the centre, and Ced 111 near IRN. The stellar clusters near 
Ced 112 and Ced 110 are associated with massive molecular 
clumps. Besides the immediate surroundings of Ced 112 and 
Ced 110, also the most prominent molecular filaments in the 
northern and central parts of the cloud, i.e. the Northern arm, 
the 40' Arc, and the central filament near Ced 110, are lined 



with young stars. This suggests that these stellar groupings and 
large scale molecular structures have a common origin. 

In contrast to the situation near the two other reflection neb- 
ulae, clumps around Ced 111 and the IRN belong to the least 
massive in the cloud, and it seems that most of the molecular 
material around young stars has dispersed on this side of the 
cloud. 

In the southern parts, south of the galactic latitude — 16?05, 
the surface density of stars is clearly lower than elsewhere 
in the cloud. In particular, there is only one ISOCAM mid- 
IR excess member candidate source and significantly fewer 
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Fig. 6. Clumps in the central and northern part of Cha I. 
Output of the cLsurf routine. Each clump identified with the 
Clumpfind programme is shown as an isointensity surface in 
the (x, y, v) space. Different colours are used to help in distin- 
guishing between separate clumps. The intensity level chosen 
for this figure is T\ = 1 K. The x and y axes represent the off- 
set in arcminutes from the map centre (R.A. = ll h 09 m 13?3, 
Dec. = -77°29'12"), and the v axis represents the LSR veloc- 
ity in km s~ 1 . Note the large velocity gradient between the two 
northernmost clumps 1 and 2 at y sa 55' (here on the right, the 
viewing angle makes the y-offsets look larger) and the excep- 
tionally large velocity of clump 8 in the centre at y w 10' (here 
on the left). 

NIR member candidates than in the north. It should be noted, 
however, that the ISOCAM observations did not go below 
ru — 77°47'. The two parallel southern filaments were de- 
tected in FIR surface emission by IHaikala et al. 19981 and 
IToth et al. 200 Considering the small number of NIR and 
mid-IR sources, these structures are likely to be cold and their 
clumps are possibly in their pre-star formation stage. (The same 
might be true for the starless clump (12) south of Ced 112, even 
though it lies next to a string of young stars on its eastern side.) 

The large number of newly born stars and their concentra- 
tion around filaments seen in the molecular line map indicate 
efficient compression of gas and subsequent star formation. 
According to the hydrodynamical models (e.q. Klessen 2001 1, 
these features are characteristic of cloud fragmentation in the 
conditions of little turbulent support (decaying turbulence), or 
turbulence driven by large-scale shocks. In contrast, turbulent 
support driven on small scales leads to inefficient and dispersed 
birth of stars. The modest C 18 column densities, the fairly 
large fraction of visible stars amongst the associated YSOs, and 



the fact that only two protostar candidates have been found in 
the cloud, suggest on the other hand that we are witnessing the 
aftermath of star formation. The total mass of the young stars 
associated with the cloud has been estimated to be about 120 
M Q JMizuno et al. 19991 . which corresponds to about 50% of 
the mass of high column density gas traced by C 18 0. 

3.3. Orientation of the magnetic fields 

Near-infrared (77-band) polarimetry of heavily reddened 
background stars and embedded objects was used by 
[McGregor et al. 1994| to study the magnetic field structure in 
Cha I. The stellar polarization due to aligned dust particles in 
the direction of the cloud is illustrated in Fig. [8] The data is 
from |McGregor et al. 199^ and lWhittet et al. 19941 The polar- 
izations of some of the Cha I members are influenced by in- 
trinsic polarization due to scattering in the circumstellar dust 
or the surface of the molecular cloud. Depending on the ge- 
ometry, this can lead to depolarization and/or rotation of the 
polarization direction. A notable example is the IRN (polariza- 
tion not plotted in Fig. [8jl for which the H band polarization 
is 33% and the polarization angle differs by tens of degrees 
from the cloud average ( McGregor et al. 1994). The small fore- 
ground reddening in the direction of Cha I and the high galac- 
tic latitude makes it certain that the observed high polariza- 
tions are in major part due to dust in the Chamaeleon region. 
These observations probe magnetic fields in the obscured parts 
of the cloud and have thus a direct bearing on structures re- 
vealed by the present C 18 mapping. The polarization field 
is well ordered throughout the cloud, changing its orientation 
from nearly E-W in the northern part to NW-SE in the centre 
and in the south. This orientation, which is consistent with the 
large-scale polarization field in the region JWhittet et al. 1994i . 
is roughly perpendicular to the cloud axis, the dense filaments 
seen in C 18 and presumably also to the magnetic field in the 
region. As pointed out by |McGregor et al. 1994| the assumed 
orientation of the magnetic field suggests that it has steered 
the cloud collapse preferentially along the field lines. That also 
dense filaments seem to be oriented at right angles to this mag- 
netic field lends support to this view. 

It should be noted, furthermore, that turbulent fragmenta- 
tion models predict that the resulting magnetic field in dense 
filaments or sheets is parallel to the direction of elonga- 
tion, because this magnetic field component is amplified in 
shocks ( Pad oan & Nord lund 2002 ). The iJ-polarization per- 
centage versus visual extinction in Cha I seems to saturate at 
Ay ~ 7 (see Fig. 4 of McGrego r et al. 1994) , suggesting that 
the grain alignment weakens in the densest parts (see the dis- 
cussion in Good man et al. 1992b . 

4. Comparison with NANTEN cores 

The integrated C ls O( J =1 — 0) intensity map presented in 
Fig.^shows the same general outline as the one taken with the 
4-m NANTEN telescope (Mizu no et al. 19 99 1. As expected, 
the higher spatial resolution of the SEST allows one to re- 
solve more details. The cores detected by Miz uno etal. 1 999 
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are all seen in the SEST data but these cores are fragmented 
into smaller entities, both spatially and in velocity. 

The correspondence between 'cores' identified by 
Miz unoetal. 1999l in the NANTEN map, and 'clumps' identi- 
fied with the Clumpfind programme in the SEST map is most 
clear in the northern arm or Cha I North. Our clumps 1 and 2 
in Cha I North with a total mass of about 10 M Q correspond to 
the centre of NANTEN core 3 (N3), for which they derived the 
mass 29 Mq. The velocity difference between clumps 1 and 
2 is clearly visible in Fig.|6](at y w 55')- Clump 12 (12 M ) 
about 15' south of Cha I North corresponds to core N4 
(16 Mq). The SEST map shows, however, several subsidiary 
peaks also in the north, and especially the windswept shape of 
the northern filament is clearly visible there. 

In Cha I Centre and in the southern part of the cloud, where 
the structure is more complex, the two NANTEN cores (6 and 
5) decompose into several clumps. Towards the absolute inten- 
sity maximum of the cloud, the peak of core N6, two clumps 
lie in the line of sight at different radial velocities (clumps 8 
and 17 at y 10' in Fig. [6). The existence of two velocity 
components was already noted bv lMizuno et al. 19991 Besides 
these two, core N6 comprises two major clumps (3, 19) and 
has several smaller clumps on its western side. The total mass 
of the four most massive clumps is 30 Mq, which roughly cor- 
responds to the mass derived for Mizuno's core 6, 39 Mq. 

The centre of core N5 in the south correspond to our clumps 
5 and 6 which belong to one of the two parallel filaments 
in the south . The procedure we have used identifies half a 
dozen of clumps in this filament with a total mass of 30 Mq, 
which is well below the mass derived for core N5 (82 Mq) 
by Mizu no~etal. 19991 The most massive individual clumps in 
the two southern filaments are numbers 6 and 16, respectively, 
which also are among the most massive in the whole cloud. 

5. Notes on individual regions 

5. 1. Cha i North 

The Clumpfind program finds in this region three major 
clumps, 1, 2, 10. Clump 1 is in the direction of the B9V 
star HD97300, which is the illuminating source of Ced 
112. Clump 2 corresponds to the opaque core described in 
IJones et al. 19851 It contains the Class protostar candidate 
Cha-MMS2 (Rei purth et al. 1996) , which is a possible driv- 
ing source of the molecular outflow detected in this re- 
gion JMattila et al. 1989i . Cha-MMS2 lies in the apex of the 
blueshifted outflow lobe. Numerous mid-IR sources with mid- 
IR excess are seen projected in the line of sight to clumps 1 and 
2 ( IPersi etal. 19991 IPersi et al. 20001 . One of these sources, 
ISOCAM-ChaINa2 has a SED characteristic of a Class I source 
and has been proposed by Persi et al. 19991 to be the possible 
exciting source of the molecular outflow. 

Clump 12 (the centre of N4), which lies 15' south of clumps 
1 and 2 is devoid of mid-IR sources and does not seem to 

1 Note: The coordinates of clump 5 in Table 1. of 
Mizu no et al. 19991 are not correct. The correct coordinates are 
I, b = 297?23, -16? 13 (Mizuno 2003, private communication.) 



be connected with star forming activity. This core was de- 
tected also in the ISOPHOT 170 /jm serendipity survey of 
the Chamaeleon clouds JToth et al. 200 1, which indicates the 
presence of very cold dust in this clump. 

5.2. Cha i Centre 

The prominent intensity maximum near Ced 110 (N6) is split 
into four major clumps: 3, 8, 17 and 19. Projected on the sky, 
clump 8 partially overlaps with clump 17, but it has a larger 
radial velocity (see Fig |6). This region contains a cluster of 
low-mass young stellar objects which have been studied exten- 
sively in recent years (e.g. lPrusti et al. 1991IIPersi et al. 2 000; 
IPersi et al. 200T1 ILehtinen etal. 20011 ILehtinen et al. 20031 
The dense dust ridge detected at 200 with ISOPHOT 
(ILehtinen et al. 200TT > coincides with the central parts of 
clumps 3 and 8. The Class candidate Cha-MMSl 
( Reipurth et al. 1996) is probably embedded in clump 3, 
close to its north-eastern boundary. |Reipurth et al. 19 96 sug- 
gested that Cha-MMSl is the central source of the bipolar 
molecular outflow discovered by Ma ttila et al. 19891 whereas 
ILehtinen et al. 2003l regarded the Class I infrared source IRS 4 
as a more likely candidate. In the latter case the outflow would 
originate between clumps 3 and 17. 

The clumps 60, 35, 19, 33, 18, 27 and 25, together with 
the clumps west of the Ced 111 region, 22, 23 and 30, form a 
nearly continuous arc of 40'. The clumps in the arc are however 
readily separated both spatially and in velocity (Fig. [2). 

5.3. Ced 1 1 1 region 

Besides the Ced 111 reflection nebula, the southeastern edge 
of the cloud contains the IRN and several pre-main-sequence 
stars (eg. Sch wartz 199U . Only four clumps, 22, 23, 30 and 38 
are located in the region which has formerly been a centre of 
active star formation. It seems that this process has consumed 
and dispersed most of dense material of the ambient cloud. The 
clump 38, which is detected at a low level, is associated with 
the IRN. 

5.4. Cha i West and South 

The clumps 4, 7, 9, 13, 14 and 15 to the West of Ced 111 
form a massive condensation of clumps. This conglomera- 
tion of clumps is probably associated with the Very cold 
Core 5 (VCC5) in IToth et al. 20001 The two parallel elon- 
gated structures to the the South are clearly seen in the far- 
IR ( IHaikala et al. 19981 ITothet al. 2000> . IT6th et al. 2000l des- 
ignated these structures as VCC4 and VCC3. Similar to VCC4 
also the two other separate into smaller units. The clumps 5 
and 6 correspond to VCC4 and clumps 16 and 24 to VCC3. 
None of these structures contain IRAS point sources with 
colours typical of newly born stars, and seem to have no ac- 
tive star formation. The ISOCAM mid-IR Chamaeleon map- 
ping (IPersi et al. 20 00 1 covered only the northern part of these 
regions. 
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6. Distribution of clump masses and their 
stability 

The distributions of the derived clump masses, the ratios of 
the potential to kinetic energies, and the ratios of the thermal 
to turbulent kinetic energies are shown in Fig. [9] The clump 
masses range from 0.6 to 12 M Q with the median 2.7 M . The 
l-Epotl/^kin ratios are between 0.1 and 1.2, and the median 
ratio is 0.4. The ratio is close to unity for the most massive 
clumps. For most clumps the turbulent energy dominates over 
the thermal energy. 

As the condition of the virial equilibrium (in the absence of 
an external pressure) is \E pot \ / = 2, the result indicates 
that, strictly speaking, none of the clumps identified here is a 
gravitationally bound entity. They are either dissolving or sta- 
bilized by the outside pressure due to interclump turbulence or 
by magnetic fields. 

We may estimate for each clump a hypothetical external 
pressure, P e xt> required to balance the 'extra' turbulent motions 
using the virial equilibrium equation: 

2E kin + E pot = 4nR 3 P cxt . (1) 

The distribution of these 'equilibrium pressures' is shown 
in Fig. [H)] The distribution is asymmetric, peaking towards a 
minimum which lies close to the value 10~ 12 Pa (P/k ~ 7 10 4 
cm _3 K). From observational statistics it is likely that the exter- 
nal pressure may be sufficient to stabilize some of the clumps, 
but not to force any of them to collapse. In case the pressure 
of the interclump medium is roughly constant across the cloud, 
the minimum of the derived P cxt values should correspond to 
the actual interclump pressure. One can see in Fig. ^| that in 
fact for about 40 clumps P cxt is close to 10~ 12 Pa. We assume 
therefore that this value represents the typical interclump pres- 
sure, and that those 40 clumps are in the virial equilibrium. The 
mass distribution of these 40 clumps is presented as hatched 
columns in Fig. [9] It is similar to that of the original sample, 
but these clumps have less marked turbulent motions (see the 
bottom panel of Fig. [9jl. In case the clumps are stabilized by 
the magnetic field, the field strength is 16 fiG, as can be ob- 
tained from the following formula for the magnetic pressure: 
P m = 5 2 [AiG]/(87r)10- 13 Pa. 

The clump mass spectrum can be, in principle, 
used for studying the mechanisms of cloud frag- 
mentation, and, if sufficiently small scale structures 
can be resolved, also for predicting the stellar initial 
mass function, IMF (see e.g. Padoan & Nordlund 2002 
Ball esteros-Paredes & Mac Low 20021 IKlessen 20011 . The 
spatial resolution and the RMS noise level of the present 
data are clearly not sufficient for studying the complete mass 
spectrum down to clumps at < O.5M . Stars may form in the 
interiors of the larger, quasi-stable clumps identified here, but 
the relation between their masses and the stellar mass function 
is beyond the scope of this paper. 

We discuss briefly, however, the possibility to understand 
the mass spectrum of the hypothetically virialized clumps 
in terms of the turbulent fragmentation model presented in 
IPadoanetal. 19971 and iPadoan & Nordlund 20021 These stud- 
ies predict 1) the mass distribution of dense cores, 2) the proba- 



bility density function of the gas density, and 3) the subsequent 
mass distribution of collapsing clumps arising from supersonic 
turbulence. 

In the model of Pa doan & Nordlund 20 02 the mass distri- 
bution of all clumps follows a power law depending on the 
power spectrum of the turbulence: N(m) oc to~ 3 / ( 4 ~W, where 
j3 is the spectral index of the turbulence (see their Eqs. (5) and 
(18)). Using their assumption that the density and mass distri- 
butions are statistically independent, one can derive the distri- 
butions of the Jeans' masses, <1>j (m), or alternatively, the distri- 
bution of 'equilibrium masses' <£>e(w), under the condition of 
constant external pressure expressed in Eq.^ in the latter case 
the equivalent of Eq. (24) of Padoan & Nordlund 2002 can be 
written as: 

$ E (m)am _3/(M) / g(m)dm, (2) 

where g(rh) is the probability density function of equilibrium 
masses, and m is the variable of integration. The function g is 
derived below. The observed distribution of pressure balanced 
clumps (in the top of Fig. [9} should be a sample of the distribu- 
tion $e- 

The /3-index should be reflected in the observed size- 
linewidth relation, Av oc R a , since a — (/3 — l)/2 accord- 
ing to Eq.(13) of Pad oan & No rdlund 2002; The expected val- 
ues range from 1 .6 (for incompressible turbulence) to 2 (shock 
dominated turbulence; Larso n 1979IILarson 198U . The clumps 
identified in the present study show a very weak dependence 
between the line width and size, and a large scatter. A least- 
squares fit to the identified clumps gives a = 0.16±0.09 (in- 
stead the usual a w 0.5) with a small correlation coefficient 
of r — 0.2. The AOS channel width and the mapping step size 
correspond to 0.12 kms -1 and 0.044 pc, respectively. These 
values correspond to a large fraction of the total range of the 
observed line widths (from 0.4 to 0.9 kms" 1 ) and sizes (from 
0.08 to 0.22 pc). Therefore, even if a correlation between the 
line widths and clump sizes exists in the cloud, the present data 
are not suitable to make a meaningful fit when the errors are 
taken into account. The qualitative results presented here do 
not depend, however, strongly on the actual value of (3. Since its 
determination from observations is furthermore subject to large 
uncertainties, we assume in the following the value 2 valid for 
compressible turbulence. 

We next derive the function g(m) appearing on the right 
hand side of Eq. |2] According to the turbulent fragmentation 
models of Padoan et al., the probability density distribution of 
the local density parameter, p(x), where x = n/n, and n is the 
average number density in the cloud, follows approximately 
a lognormal function. The probability density distribution per 
unit mass, i.e. the mass function at any given density, f(x) = 
xp(x), can then be written as 




(see Eqs. (1) and (7) of Padoan ^t al. 19971 . As pointed out 
in the latter study, the properties of the lognormal function 
and the selection of the variable x such that x — 1 imply 



L. K. Haikala et al.: C ls O (J = 1 - 0) in Chamaeleon I 



11 



that the mean, Inx, and the standard deviation, a\ nx , are re- 
lated by In a; = —Ofa x /2- The distribution is thus charac- 
terized by a single parameter, a\ nx . This in turn depends on 
the rms Mach number, A4, of the turbulent flow according to 
a? nx « ln(l + 0.257W 2 ) JPadoanetal. 19971 Eq. (9)), i.e. in- 
homogeneity grows with the speed of the flow. 

Turning to a small element of the cloud, we examine an 
isothermal, homogeneous clump formed within it. From Eq. \l\ 
one can derive the following expression for the mass of a spher- 
ical, virialized clump as a function of the density parameter x: 



m(x)=ax- 1/2 (b-cx- 1 ) 3 / 2 



(4) 



where 



5_ 

G 

kT 

1 

l L 

Pcxt 



.3/2 



( 3 



\4.7r fin 



1/2 



(5) 



where G is the gravitational constant, k is the Boltzmann con- 
stant, T is the gas kinetic temperature, fi is the mean molec- 
ular mass, and a v is the turbulent velocity dispersion inside 
the clump. For Jeans' masses the corresponding function can 
be written as mj(x) = (|) 3 / 2 a6 3 / 2 x -1 / 2 with the same def- 
initions of a and b as in Eq. [5] One can see from Eq. |4] that 
a minimum density, xo = c/b, is required for an equilibrium 
mass. This condition is equivalent with the requirement that 
the internal pressure, P; nt = nfi(kT//j, + a 2 ), must be equal to 
or larger than P C xt- The equilibrium mass reaches a maximum 
at x = 4a; = 4c/ 6, where m = (3 3 / 2 /16) a b 2 c" 1 / 2 , and 
decreases monotonically thereafter. At the maximum m(x) = 
mj(x) and E pot — — f-Ekin- This turning point signifies the 
density at which gravity starts to dominate over external pres- 
sure as the counterforce against the internal pressure. 

The functions m(x) and m,j (x) are plotted in Fig. II lh. The 
kinetic temperature is assumed to be 10 K, and the average gas 
density and velocity dispersion derived for the C 18 clumps 
(n — 5500 cm -3 , a v — 260 ms"l) are used to estimate the 
values of the parameters a, b and c. It can be seen from this 
figure that m{x) exceeds mAx) when x > Ac/b. According 
to the definition of Jeans' mass this means that the equilibrium 
masses on the right hand side of the maximum will collapse if 
the external pressure increases. 

The probability density function of equilibrium masses, 
g{m), can be derived from Eqs. [5] and 0] For each mass be- 
low the maximum derived above there are two possible values 
of the density parameter x, say x\ (low density, < Axq) and X2 
(high density, > 4xo, 4xq corresponds to about 10 4 cm" 3 with 
the parameters used for Fig. II 1>. 

If we denote by M and X the random variables describing 
the mass and the density parameter, respectively, the probabil- 
ity distribution function of the mass can be expressed with the 
aid of the density as follows: 



P(M < to) = P(X <xi) + P(X > x 2 ) 



f{x)dx 



lx 2 (m) 



f(x)dx 



(6) 



The probability density function, g(m) can then be obtained 
from the latter by derivation: 



flH = /Oi(m)) 



dxi (to) 



dm 

f{x 1 {m)) f{x 2 (m)) 



f(x 2 (m)) 



dx 2 (to) 
dm 



m'{x\{m)) m'(x2[m)) ' ^ 

where function / is defined by Eq.|3j x\ (m) and x 2 (to) are the 
inverse functions of m(x) in the regions xq < x < Axq (low 
densities) and x > 4xq (high densities), respectively, and to' 
is the derivative of function m. The expressions of x\ (to) and 
x 2 (m) that follow directly from Eq. 0]are rather intricate and 
they are not shown here. The functions X\{m) and x 2 (m) are, 
however, plotted in Fig. II lb. together with the inverse of Jeans' 
mass, Xj(m). 

Estimates for the parameters a, b, c and (3 in Eqs. |4]and|2] 
can be derived from observations or, as in the case of (3, from 
theoretical predictions. The mass distributions of equilibrium 
and collapsing clumps predicted by the model of Padoan et 
al. depend then ultimately on the gas density distribution de- 
scribed by the parameter a\ n x - It turns out that a small disper- 
sion in densities (cri n x < 1-2) makes the Jeans' mass distribu- 
tion peak towards high masses, M > 10M Q , which contradicts 
the observed stellar mass distribution in Cha I and in other dark 
clouds. A situation where both Jeans' masses and equilibrium 
clump masses peak below 1 M require a wide range of densi- 
ties, and consequently, a large value of <j\ n x , typically in excess 
of 2.5. These high values in turn imply large turbulent Mach 
numbers (A4 > 45, corresponding to 8.6 kms" 1 at 10 K or 30 
kms^at 100 K: IPadoan et al. 19971 Eqs. 4 and 5). This situ- 
ation may have resulted from the passage of a powerful shock 
wave, e.g. one caused by a collision with an expanding bubble 
driven by supernova explosions and stellar winds from an OB 
association. 

The density distribution for a slightly smaller density dis- 
persion, <j\ nx = 2.2, and the resulting mass distribution are 
illustrated in Figs. II It and d, respectively. The latter Figure 
shows the probability density function of 'equilibrium masses', 
$e(w), obtained from Eqs.|2]and0 and the corresponding dis- 
tribution of Jeans' masses, $ j (to),w together with the 'mother' 
distribution of all clumps, assumed to follow the power law 
to" 15 . The value of a\ nx is chosen so that <£>e(ti) decreases 
towards higher masses, but the function still shows the charac- 
teristic peaking towards zero, which cannot be distinguished 
at larger values of a\ nx . The peak near zero is the contri- 
bution of low-density, pressure balanced clumps, whereas the 
smooth curve peaking near 3M© represents dense, gravitation- 
ally bound clumps. The latter distribution resembles that of the 
Jeans' masses, but is flatter. 

The clumps detected in this study are likely to bath in 'inter- 
clump' gas with the density just below the critical density of the 
C ls O( J =1 — 0) line, i.e. about 10 3 cm -3 . Assuming that the 
external pressure 10" 12 Pa is caused by turbulence in this inter- 
clump gas, and that the gas temperature remains at 10 K also 
there, we find that the required turbulent velocity dispersion 
is 0.5 kms" 1 or that the rms Mach number is about 3, i.e. 15 
times lower than the one consistent with the observed mass dis- 
tribution according to the model of Pado an & Nord lund 2002 
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Fig. 11. a) The 'equilibrium mass' m(x), i.e. the clump mass as a function of the density parameter x = n/n (solid line). 
The physical parameters used here correspond to values estimated for the C 18 clumps in Cha I: P cxt = 1 1CP 12 Pa, n = 
5.5 10 3 cm" 3 , T = 10 K, a v = 260 ms -1 . The Jeans' mass, mj, for these conditions is plotted with a dashed line, b) Density 
parameter a; as a function of mass, i.e. the inverse of m(x) in two regions x±(m) (dotted line, low densities) and x% (m) (solid 
line,high densities). The corresponding function for Jeans' masses, Xj, is plotted with the dashed line, c) The assumed lognormal 
density distribution, p(x), with the standard deviation a\ nx — 2.2. d) The probability density function of equilibrium masses, 
$e(w)- The distribution of the Jeans' masses, <&j(to), is plotted as a dashed curve. The dotted curve represents the assumed mass 
distribution of all dense cores (not in scale). It follows the power law m 5 , which is consistent with compressing turbulence. 



The discrepancy may be partly explained by the hierarchi- 
cal structure of turbulence. If the rms Mach number 3 is char- 
acteristic for the length scale 0. 1 pc (corresponding to the typ- 
ical size of a clump), M. = 45 (a v w 8.5 kms" 1 ) should 
then characterize turbulence in the length scale 20 pc, assum- 
ing that the scaling law is cr v oc L 0,5 valid for compressible tur- 
bulence. This length scale corresponds to the size of the whole 
Chamaeleon dark cloud complex including the clouds I, II and 
III. The time scale, L/a, associated with it is 3 Myr, meaning 
that in this time the large-scale turbulence should have decayed 
into motions on small scales. The ages of the star associations 
do not contradict with this time scale. As discussed in Sect. 
3.2., a large fraction of the stars associated with the Cha I are 
visible pre-main-sequence stars, which indicates that at least 
a time on the order of 10 6 years has passed since the initia- 
tion of star formation and the possible violent event triggering 
it. On the other hand, if the quoted scaling law is valid, one 
could expect to see a larger dispersion in the radial velocities 
of the clumps across the cloud, i.e. over a projected distance 



of about 4 pc, than the value 0.35 km s 1 mentioned earlier in 
this Section. 



The exercise performed above shows that the origin of the 
observed clump mass distribution may be explained by the tur- 
bulent fragmentation model of Pa doan & Nordlund 20 02 if tur- 
bulence has initially been hypersonic, but cascaded down to 
much smaller scales and speeds by now. The deduced present- 
day velocity dispersion in the interclump gas suggest that the 
length scale has initially been larger than the whole Cha I 
cloud. According to numerical simulations turbulence acting 
on large scales gives rise to extensive filaments of dense gas 
(e.g. Klessen 2003 1, and the structure observed in Cha I could 
be readily understood. However, as discussed in Sect. 3.3 the 
well-ordered magnetic field is not consistent with this picture. 
Moreover, it is not clear from simulations if the velocity disper- 
sion across the cloud can become as low as observed in Cha I 
in the turbulent fragmentation process. 
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In view of the fact that the region studied here may cover 
only a small fraction of the structure formed in this process, it 
is questionable how well it can reflect the global properties of 
the flow. We therefore do not proceed to a quantitative anal- 
ysis, but intend to present such in a subsequent paper with a 
more varied sample of clumps, obtained by combining results 
from different clouds belonging to the same complex and using 
different clump identification algorithms. 

7. Conclusions 

The C 18 (J = 1 - 0) map of the Chamaeleon I cloud 
with an angular resolution of ~ 1' shows a network of fila- 
ments with the likeness of some recent simulations of molec- 
ular cloud fragmentation driven by large scale turbulence (e.g. 
IKlessen 20 01 1. These filaments seem to have, however, pref- 
erential directions, which change from N-S in the north to 
NE-SW in the other parts. 'Cores' identified in the previ- 
ous NANTEN map ( Mi zunoetal. 19991 break up into a num- 
ber of substructures, which we call 'clumps'. We have ana- 
lyzed the small scale structure using the 'Clumpfind' routine 
( Williams et al. 1994 1 and identified about 70 separate clumps 
in the (x, y, v) data cube. The RMS noise level of these data 
limits the detectable clump spectrum to masses above 0.5M Q . 
The inspection of the clump masses and kinematics shows that 
none of these clumps is a gravitationally bound entity, unless 
a large fraction of the mass is hidden from C 18 observations 
because of depletion. However, a low external pressure of about 
10~ 12 Pa (P/k ~ 710 4 cm" 3 K) would be sufficient to bring 
some 40 clumps to a balance. This kind of pressure could be 
caused by turbulent interclump medium (n ~ 10 3 cm~ 3 ) with 
a velocity dispersion of a v = 0.5 kms -1 or a magnetic field 
with a strength of 16 /iG. 

Using the analytical model for turbulent fragmentation pre- 
sented in Pad oan et al. 19971 and IPadoan & NoFd lund 2002|we 
derived a theoretical mass spectrum for clumps in virial equi- 
librium by help of external pressure. It turned out that large 
turbulent Mach numbers are required to produce a mass spec- 
trum peaking at low masses as observed, and that correspond- 
ing turbulent velocity dispersion clearly exceeds the interclump 
velocity dispersion quoted above. A possible reason for the dis- 
crepancy is that the turbulence in the cloud has decayed since 
the hypothetical equilibrium clumps formed. Alternatively all 
condensations seen in the C 18 map are just transient struc- 
tures. In any case, gravitationally bound protostellar cores, such 
as Cha-MMS 1 and Cha-MMS2 probably are, represent subcon- 
densations of C 18 'clumps'. 

Although the clump mass distribution may be understood in 
terms of turbulent fragmentation models, one observational re- 
sult is not easily explainable by them, namely, the well-ordered 
magnetic field structure in the cloud. The magnetic field di- 
rection as determined from the polarization position angles of 
highly reddened background stars is perpendicular to the gen- 
eral orientation of dense filaments, suggesting a magnetically 
controlled collapse (McGregor et al. 1994). 

We consider it possible that different processes are taking 
effect on different scales. Although large scale magnetic fields 
are likely to have influenced the formation of the filaments, 



their fragmentation into smaller clumps may have been dom- 
inated by chaotic motions. In the sequel, we intend to perform 
direct comparisons between the observational data (spectra and 
column density distributions) and simulation results in order to 
determine the characteristics of the turbulence and thereby gain 
better understanding of its role in the cloud fragmentation. 
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Fig. 7. Cha I association members and candidate mem- 
bers projected on the integrated C ls O intensity map of 
the cloud. Sources of the coordinates are ISchwartz 19911 
(Optical candidate members, black triangles), P ersi et al. 2 000 
(known members with NIR excess, red squares, mem- 
ber candidates with NIR excess, open red triangles), 
|Cambresy et al. 1998| (new Young Ste llar Object candidates, 
asterisks) and |G6m ez & Kenyon 2001 (candidate pre-main- 
sequence stars, open circles). The locations of the reflection 
nebulae Ced 112, Ced 110 and Ced 111 (from north to south) 
are indicated with open squares. The locations of the Class 
protostar candidates Cha-MMS 1 and Cha-MMS2 are indicated 
with filled circles. 
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Fig. 8. Polarization vectors of field stars and Cha I members 
( Wh ittet et al. 19941 |KlcGregor et al. 1994 1 superposed on the 
integrated C ls O intensity map of Cha I. 
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Fig. 9. Top: Distribution of the clump masses in Cha I. Middle: 

The ratios of the gravitational potential to the kinetic energies, 
l-EpotlZ-Ekinl- The kinetic energy estimates contain both the 
internal thermal energy and the energy of the turbulent mo- 
tions. Bottom: The ratios of thermal to turbulent energies. The 
hatched columns correspond to clumps which possibly are in 
virial equilibrium (see text). 
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Fig. 10. The distribution of external pressures, P ext needed to 
stabilize the identified clumps. 



